How to Measure Forces when the Atomic Force Microscope shows Non-Linear 

Compliance 
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A spreadsheet algorithm is given for the atomic force microscope that accounts for non-linear 
behavior in the deflection of the cantilever and in the photo-diode response. In addition, the data 
analysis algorithm takes into account cantilever tilt, friction in contact, and base-line artifacts such 
as drift, virtual deflection, and non-zero force. These are important for accurate force measurement 
and also for calibration of the cantilever spring constant. The zero of separation is determined 
automatically, avoiding human intervention or bias. The method is illustrated by analyzing measured 
data for the silica-silica drainage force and slip length. 



I. INTRODUCTION 

Scanning probe microscopy has revolutionised surface 
science by enabling the production of images of surfaces 
with molecular resolution. In one technique a topo- 
graphic map is produced based on the height adjustment 
of a piezo-drive required to maintain a constant deflection 
of the cantilever probe during a raster scan. The original 
atomic force microscope used tunneling currents to detect 
the deflection of the cantilever^ This was soon modified 
to use a light lever to detect the deflection, which had the 
advantage of not requiring vacuum conditions. 2 '- It also 
allowed quite large deflections to be measured, and differ- 
ent imaging techniques to be developed, such as constant 
height imaging. For such techniques to be quantitative, 
the light lever signal (in volts) has to be converted into 
the deflection of the cantilever (in nanometers). 

This issue of quantitatively calibrating the light lever 
received added impetus with the further modification of 
the atomic force microscope in what has come to be called 
colloid probe force microscopy^ In this technique a col- 
loid sphere of measured radius (R s» 10-20 /xm) is glued 
to the end of the cantilever spring instead of the sharp 
tip used for imaging. The object is to measure the so- 
called surface force between the substrate and the probe 
as a function of separation. The surface force is just the 
spring constant times the cantilever deflection, and so if 
the light lever is properly calibrated, the measurement 
may be performed with molecular resolution. 

The light lever itself is usually made from an optical 
beam reflecting off the back of the cantilever spring onto 
a split photo-diode. A change in angle of the cantilever, 
due, for example, to a change in the force on the probe, 
causes the light beam to move across the face of the 
photo-diode. The consequent change in voltage differ- 
ence between the two halves is measured and taken to be 
proportional to the change in cantilever angle. By using 
the piezo-drive to press the cantilever against the hard 
substrate, the proportionality constant is obtained as the 
slope of the voltage versus distance signal. 

To make this clear mathematically, let the measured 
constant compliance slope be 
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FIG. 1: The raw photo-diode voltage versus the piezo-drive 
displacement. The solid curve is measured extension data, 
the dashed line is the tangent to the contact region at first 
contact, and the dotted line uses the average slope in con- 
tact. The inset shows the analysed force versus separation, 
with the solid curve resulting from the non-linear analysis 
(see text) and the dashed and dotted curves resulting from 
the conventional linear analysis using the first contact slope 
and the average slope, respectively. The source of the mea- 
sured data is Ref. [fj; a summary of the experimental details 
is given in Ref. 0. 



Here V is the photo-diode voltage and z p is the piezo- 
drive position. Letting Vb be the base-line voltage far 
from contact, the surface force is 



F(z p ) = -ko/3- 1 {V(z p )-V h } 1 



(2) 



AV 



(1) 



where fco is the cantilever spring constant. The force 
goes to zero at large separations, which is the base-line 
region. This result invokes the fact that in hard con- 
tact the change in tip position is equal and opposite to 
the change in piezo-drive position,— Az t = — Az p , and 
the force on the cantilever is F — koz t . (The simple 
presentation given here ignores a number of important 
linear effects such as base-line drift, friction in contact, 
cantilever tilt, and non-negligible base-line force. These 
and non-linear effects will be included in the more so- 
phisticated linear and non-linear analysis below.) The 
separation is 

h = z p + zt + const 
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= z p + k 1 F cxt (z p ) + const, (3) 

where the constant is chosen so that ft, = in contact. 

It ought to be clear from the above how important the 
calibration factor /3, which is the constant compliance 
slope, is to the quantitative measurement of surface forces 
with the atomic force microscope. However it is not un- 
usual for the photo-diode voltage versus piezo-drive po- 
sition curve to be non-linear in the contact region. A 
typical example is shown in Fig. [TJ It is emphasised that 
the data in the figure were obtained for hard surfaces and 
so the curvature evident is not due to elastic deformation 
of the probe or substrate. A curved constant compliance 
region such as that in Fig. [TJ creates several problems. 
At a minimum, the calibration factor /3 is not unique; 
the slope depends upon where in the contact region it 
is measured, which introduces quantitative uncertainty 
into the linear analysis, (compare, for example, the re- 
sult that uses the slope at first contact (dashed line and 
curve) with the result that uses the average slope (dot- 
ted line and curve)). Ambiguity also arises because the 
contact behavior differs between extension and retraction 
(not shown). Worse, a non-linear contact region contra- 
dicts the fundamental assumption of a linear response 
and raises questions about the conventional linear anal- 
ysis that is used to quantify surface force measurements. 

The inset to Fig. [1] graphically illustrates the prob- 
lems with the linear analysis. It can be seen that it 
gives grossly unphysical non-zero separations in the con- 
tact region. In this case the surfaces are known to be 
rigid so the results is unambiguously unphysical. In other 
cases, where the surfaces are either of unknown rigidity 
or known to be soft, this artefact of the linear analy- 
sis would be misinterpreted as the elastic deformation 
of the material. Not only the elasticity but any useful 
physical information in the contact region is precluded 
by the linear analysis. As will be demonstrated below, 
the non-linear analysis can be used to obtain reliable val- 
ues for properties like the friction coefficient, roughness, 
and topography of the contact region. 

It is not only in contact that the linear analysis can 
fail. If a large surface force is present, then the linear 
analysis of the data introduces quantitative errors into 
the values of the surface force in the non-contact region. 
This can be a particular problem if one requires reliable 
and accurate surface force measurements, or if one seeks 
small changes in forces with control parameters, or if one 
needs to quantify second order effects. In all these cases 
the linear analysis can be unsuitable, depending upon 
the extent of the non-linearity and the magnitude of the 
forces. 

The are two possible physical origins of the non- 
linearity displayed in Fig. [TJ The first possibility is that 
the cantilever deflection becomes non-linear over the rel- 
atively large range of the contact region. By cantilever 
non-linearity is meant that the four relevant quantities 
(tip position, deflection, angle deflection, and force) are 
not linearly proportional to each other. The second pos- 
sibility is a non-linear relationship between the change in 



photo-diode voltage and the change in cantilever angle. 
Since the change in voltage difference in the split photo- 
diode depends only on that part of the optical beam cur- 
rently crossing the boundary (assuming uniform sensi- 
tivity of the photo-diode; any spatial variability in the 
sensitivity will contribute further to non-linearities) , any 
variability in the spatial intensity or width of the beam 
(e.g. circular or elliptical cross-section, Gaussean inten- 
sity distribution) will give non-linear effects. 



1. Contents 



For the case of the rectangular cantilever, a complete 
linear analysis is given in SjTTJ In ijllBl the effects of base- 
line drift, virtual deflection, cantilever tilt, friction in con- 
tact, and non-negligible base-line force are accounted for. 
These are often neglected in the conventional linear anal- 
ysis. A new result is an algorithm for determining the 
zero of separation, i jll B 51 The effective spring constant 
that must be used in the linear analysis and its relation 
to the cantilever spring constant is given in B 61 In 
i jlll Al are given the non-linear equations for a rectan- 
gular cantilever that determine the deflection, deflection 
angle, vertical position, and applied force, taking into ac- 
count tilt and friction. In ElIIIBI an algorithm suitable for 
spreadsheet use is given for the analysis of experimental 
data (i.e. the conversion from raw voltage versus piezo- 
drive position to force versus separation) in the case of 
non-linear behavior. The non-linearities can arise either 
from the non-linear cantilever deflection or the non-linear 
photo-diode response, or both. For the non-linear photo- 
diode case, the non-linearity is characterized by the mea- 
surement itself and it is not necessary to know the details 
of the source of the non-linearity. This is fortunate be- 
cause unlike rectangular cantilevers, these vary between 
different models of the atomic force microscope. In £|III CI 
the case of a linear photo-diode and non-linear cantilever 
is explored numerically for the case shown in Fig. [TJ and 
it is concluded that the non-linear cantilever deflection is 
insufficient to account for the measured non-linear effects. 
In i]III D[ summarized are the equations for a linear can- 
tilever and a non-linear photo-diode, which are somewhat 
simpler than the dual non-linear case. In ijlVl these are 
applied to measured atomic force microscope data for the 
drainage force at several drive velocities. Results for the 
slip length and the drainage adhesion are obtained. The 
quantitative and qualitative differences between the lin- 
ear and the non-linear analysis of the experimental data 
are shown. 
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II. LINEAR ANALYSIS 
A. Horizontal Cantilever, No Friction 

1. Deflection 

The bending of a cantilever beam under the influence 
of fores and torques is one of the classic problems of the 
theory of elasticity. A beam of length L , and with a 
deflection x and angular deflection 9 has stored elastic 
energy^ 



U(x, 6) = ^ [3x 2 



3xL e + L 2 e 2 ] 



(4) 



The elastic parameter B = EI depends upon Youngs 
modulus and the geometric second moment of the beam; 
it will be related to the spring constant of the cantilever 
beam below. Differentiating with respect to the deflec- 
tion gives the force exerted on the end of the beam, 



F 



dU(x,9) 2B 



dx 



LI 



[6x — 3Lq6 



(5) 



and differentiating with respect to angle gives the torque 
_ dU(x,6) 



39 



^ [-3L x + 2L\i 



(G) 



These assume that the beam is in equilibrium with the 
applied forces and torques. 

Inverting these equations gives the standard expres- 
sions for the deflection and the angle in terms of the 
applied force and torque^ 



I 

2B 



2 . 
-LlF 

3 



and 



= 2B [ L o F + 2L or] 



(7) 



(8) 



2. Spring Constant 

One has to be a little cautious about assigning a spring 
constant. The above equations refer to a free, horizontal 
cantilever beam, and now the spring constant for such 
a beam with normal force and zero torque (free deflec- 
tion) will be given. It is emphasized that this cannot be 
applied to the atomic force microscope without modifi- 
cation because in that case the cantilever beam is tilted, 
the force is not entirely normal to the beam, and there 
are non-zero torques due to this and due to friction. This 
case will be handled shortly. 

For the horizontal beam with zero torque, r = 0, the 
deflection is linearly proportional to the applied force, 
x = [Lq/3B]F, from which one can identify the cantilever 
spring constant as 



k = W/L 3 . 



(9) 



In this paper the elastic parameter B will be used, as it 
is an intrinsic property of the cantilever. Usually (but 
not always; see ^11 B 61 below) any quoted or measured 
spring constant is the horizontal, free spring constant in 
the above sense, and so this equation can be used to 
convert from fco to B. 

For free deflection, t — 0, the angular deflection is 
linearly proportional to the deflection, 



2L 



x. 



(10) 



This particular proportionality constant only holds for 
the free, horizontal cantilever. For the free tilted can- 
tilever, the two remain linearly proportional to each 
other, but a different constant applies, as is derived be- 
low. 

The linear proportionality of deflection angle 9 and 
deflection x underlies the linear analysis of atomic force 
microscope. The light lever is assumed to give a change 
in voltage that is linearly proportional to the deflection 
angle, 7 = AV/A9. Since the deflection angle, deflection, 
and force, are all linearly proportional to each other by 
the above two equations, then one can get the force from 
the measured change in voltage, 



F = knx = 



2k L 



2k L 
3 7 



AV. 



(11) 



If the measured gradient of the photo-diode signal in con- 
tact is p = AV/Az p = -AV/Ax, then 



-2L n 



7 = 



(12) 



Hence F — —koP~ 1 AV. This is the conventional lin- 
ear analysis for extracting the force from of atomic force 
microscope measurements. It ought be clear it assumes 
a horizontal cantilever with no torque, neither of which 
assumption holds in practice. 



3. Spring Constant Calibration 

One of the most important issues in measuring forces 
with the atomic force microscope is the determination of 
the spring constant. This is usually the largest source 
of systematic error. The common thermal calibration 
procedure^ that is often built into the software of the 
atomic force microscope gives erroneous results, with a 
systematic overestimate of the spring constant of 15%- 
30%. The source of the error in the derivation has 
been identified and a more reliable thermal calibration 
method has been givenJ^ (The correct thermal calibra- 
tion formula for the cantilever spring constant is given 
in Eq. (|35|) below.) An even more accurate and reliable 
way of determining the spring constant is to use the long 
range hydrodynamic drainage forceJ^— In general the 
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FIG. 2: Cantilever geometry in the atomic force microscope 
(not to scale). 



with the constant Zq x * calculated to give zero separation 
at first contact, as defined below. In contact, Ah — 0, so 
that Azt = — Az p .- 

Trigonometric functions of the tilt angle occur fre- 
quently below and so it is convenient to define fixed con- 
stants 



Cq = cos 6*o, and So = sinf 



(15) 



Note that in the present geometry for the atomic force 
microscope, So ~ 0q ~ —0.2, which is small in magnitude 
and negative in sign. (All angles here and throughout are 
measured in radians; in degrees, 9q « —11°.) 



2. Friction Force 



drainage force is known exactly, at least in the large sep- 
aration regime 

Fdrain{h) = rr -, (13) 

where rj is the viscosity. It is permissible to use the piczo- 
drive velocity z p rather than the rate of change of change 
of separation in this because the deflection is small and 
its rate of change is negligible at large separations. The 
correct spring constant gives agreement between this and 
the measured force at long range. 

It should be noted that if this method of calibration is 
used in conjunction with the conventional analysis of the 
measured data (linear calibration, horizontal cantilever), 
then the spring constant that results is the effective 
spring constant fc c ft rather than the intrinsic cantilever 
spring constant ko. These are defined in the full analysis 
for the tilted cantilever with friction that is treated next. 



B. Tilted Cantilever with Friction 

1. Model 

Following earlier work ; 14 ' 15 the cantilever and probe 
in the atomic force microscope is modeled as in Fig. [2] 
The key features are the fixed tilt angle of the cantilever, 
9q < 0, such that the total angle of the cantilever is the 
sum of this and the deflection angle, 9 to t = Oq + Q, and the 
rigid lever arm, L2(#tot), which connects the cantilever 
to the point of application of the normal surface force F z 
and the lateral friction force F y , if present. The force F 
and torque r on the cantilever treated in the preceding 
section are a function of these two forces, the lever arm, 
and the tilt angle, as will now be derived. 

Note that in the figure the piezo-drive is connected to 
the base of the cantilever, so that extension corresponds 
to z p < and retraction corresponds to z p > 0.'° The 
separation between the surfaces is 

h = z p + z t + z™\ (14) 



In this work the lateral force will be taken to be due 
to friction (in contact) and it will be taken to be linearly 
proportional to the load, 

{/iF 2 , extension, contact, 
—fiF z , retraction, contact, (16) 
0, non-contact. 

This is known as Amontons law. As drawn in Fig. [21 on 
extension i p < 0, and so in contact on extension z t > 
0. This means that due to the tilt angle, y < and 
F y > 0. Since in contact F z > (at least sufficiently 
far into contact), this accounts for the sign of the first 
equality. (In general the friction coefficient is positive.) 
The opposite occurs on retraction (y > and F y < 0), 
the second equality. Out of contact there is no friction. 

The assumption that the friction force is linearly pro- 
portional to the load is a significant one. This is certainly 
the classical model of friction, at least at the macroscopic 
level. There is evidence in the atomic force microscope 
literature for— and against— such an assumption. The 
former data is perhaps the most convincing as four inde- 
pendent measurements were made (two different colloid 
probes, two different friction measurement methods). 

Immediately after the initial contact on extension, and 
at the beginning of the retraction branch in contact (the 
turn point), the probe is not moving at uniform velocity 
and so the assumed form for the friction force will pro- 
duce artifacts at these points in the analyzed data. Also, 
when the force at either first or last contact is non-zero, 
the model gives a discontinuity in the friction force and 
consequently a discontinuity in the surface force that is 
also an artefact of the simple model. 

All the following results will be given explicitly for ex- 
tension in contact. The retraction results may be ob- 
tained by the replacement fj, =>■ —fj,, and the non-contact 
results may be obtained by the replacement /i 0. The 
superscripts 'ext' (contact, fi > 0), 'ret' (contact, fi < 0), 
and 'nc' (non-contact, [i = 0) will often be used to denote 
each of the three cases. For the voltage, the piezo-drive 
position, and the surface force, which is possibly velocity 
dependent, the superscripts 'ext' and 'ret' will be used 
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also in the non-contact situation. The subscript 'c' de- 
notes a quantity in contact, and the subscript 'b' denotes 
a quantity in the base-line region far from contact. 



3. Linear Cantilever Equations 



In Eqs (|45|) and (|46|) below, non-linear expressions are 
derived for the deflection, the angular deflection, and the 
force. In the linear regime one can simply make the re- 
placement #tot Oq to obtain 



+ nS 



L Q L 2 

2B 



[S -vC ]\F z , (17) 



and 



E F, 



fiSo 



L0L2 
B 



[So - nC ] \ F, 



(18) 



Here the length of the lever arm in the undcflcctcd state is 
L 2 = L 2 (9q) = i?V2 + 2 cos 9q. In the linear regime, the 
force, deflection, and deflection angle are linearly pro- 
portional to each other. In particular it is convenient 
to define the proportionality constant between deflection 
and deflection angle from 



2L 3 [C + fiSp] + 3L 2 L 2 [S - ggbj 
3L 2 [C + fiS ] + 6L L 2 [S - (iC Q ] 
D cxt 6». 



(19) 



In general L2 <C Lq, and this and the above results could 
be expanded to linear order in L 2 / Lq. There is no great 
advantage in doing this. 

The vertical position of the tip depends upon the de- 
flection, the deflection angle, and the length of the lever 
arm, Eq. (RUl) below. The linearized form of Eq. (131))) for 
the length of the lever arm is L2(9tot) — L% — R 2 Sq8/L 2 . 
For the case of a tipped cantilever, L 2 (9 tot ) is equal to 
the length of the tip, and there is no dependence on the 
deflection angle (i.e. the term in R 2 may be set to zero). 
Linearising Eq. (|40)l below for the vertical tip position 
yields 



z t = C x + 



L 2 Sq + L 2 CqOq 



= {D ext C + L 2 So- 



+ L 2 Cq6q 



(20) 



The proportionality constant is a = dO/dzt- It has a 
different value for each of the three cases /j, > (contact, 
extension), fi < (contact, retraction), and /i = (non- 
contact). 



4- Linear Analysis of Measured Data 

The measured data in the atomic force microscope con- 
sists of the raw photo-diode voltage V(t) and the piezo- 
drive position z p (t). These are both a function of time, 
and so the voltage may equivalently be regarded as a 
function of position, V(z p ). One has in fact two sets 
of data, one for extend, V cxt (z p ), and one for retract, 
V ret (z p ). Only the equations for extension will be shown 
explicitly here. 

The tilde on the voltage is used to denote the raw mea- 
sured voltage. The raw voltage contains contributions 
from the change in angle of the cantilever and from vari- 
ous artifacts that include a constant voltage off-set, ther- 
mal drift, and virtual deflection. 12,13 Two further physi- 
cal effects have to be carefully accounted for, namely the 
drag force on the cantilever and the long range asymptote 
of the surface force. The deflection angle of the cantilever 
due to the surface forces is what is desired to extract from 
the measured voltage. The notation V(z p ) will be used 
to denote the measured voltage that has been corrected 
for these various artifacts and forces. 

In the base-line region, where the surfaces are far from 
contact, the surface force is small and in many cases neg- 
ligible. Hence almost all of the measured voltage in this 
region is due to the artifacts just mentioned. These in 
general are linear functions of position, and one can de- 
fine the measured base-line voltage on extension as 



f/cxt 



/3b xt [ 



Zpb], 



(21) 



and similarly for retraction. The base-line slope is /3§ xt = 



dV ext (z p )/dz p 



The coefficients for this are obtained 



Zpb 



by a linear fit to the measured data in an interval about 
the fixed position z p b in the base-line region. Once the 
coefficients are determined, this linear fit is applied to 
the whole measured regime, not just the base-line region 
(because the artifacts that it removes apply to the whole 
regime). With this the corrected voltage on extension is 



V cxt (z p ) = V cxt (z p ) 



V b ext (z p ). 



(22) 



This is zero in the base-line region. 

This expression removes from the raw signal not only 
the artifacts mentioned above but also the constant drag 
force. (In some cases the drag force is not constant i 12 i 13 
This effect, which can be important for cantilevers with 
a low spring constant, is not included in the present anal- 
ysis.) It also removes the linear extrapolation of the 
asymptote of surface force, 



T^CXt 



(*p) 



TTlCXt 



F 



ext f 



z pbj 



(23) 



Here the constant force is F^ xt 



derivative is F£ xt 



F cxt (h h ) and the 

= di^ ext (/i b )/d/ib, where the separa- 
tion is /ib = Zpb + ZQ Xt . (This neglects the deflection 
of the cantilever, which should be negligible; if not, add 
to the separation z t h ~ F£ xt /ko.) In almost all cases 
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this extrapolated surface force from the base-line region 
is negligible. In those cases where it isn't, it has to be 
added back, as will be done shortly. 

The contact region is where the separation between the 
surfaces is zero, and the tip moves equal and opposite to 
the piezo-drive, Az t = -Az p iS In the linear regime, the 
slope is constant and this is also called the constant com- 
pliance regime. It does not matter whether one fits the 
raw data or the corrected data because the two contact 
slopes are related by 



AV° 



AV° 



A;. 



= ft 



3cxt 

Pb • 



(24) 

If a positive voltage corresponds to a repulsive force, then 
the slope ought to be negative. 

The light lever measures the angle of the cantilever. 
The key to analyzing atomic force microscope force data 
is to calibrate the light lever by measuring the propor- 
tionality constant between angle and photo-diode volt- 
age, 



AV 
~A9' 



(25) 



This is the same on extension and retraction, and it is 
the same in contact and out of contact. This expression 
assumes a linear photo-diode, but not necessarily a linear 
cantilever. 

The value of this conversion factor follows from the 
measured slope in contact, Eq. (|24p . and the linear pro- 
portionality between angular deflection and tip position, 
Eq. (|2"0|) . Evaluating these on extension in contact one 
has 



7 



AV Azp Az t 
Az p Az t AO 



-/3 c ext K 



(26) 



One has a similar result for retraction in contact 7 rct = 
— /3£ et /a! rGt . Since this has to be a property of the light 
lever, the value of 7 cannot depend upon whether or not 
the surfaces are in contact, or whether the measurement 
is made on extension or on retraction. Hence one must 
have 7 cxt 



7 rct , or 



«(m) 



(27) 



since a oxt = a(fi) and a rot = a(— /1). The left hand side is 
a measured quantity, and the right hand side is a known 
non-linear function of fi, Eq. (|20[) . There exist sophis- 
ticated algorithms for solving such non-linear equations, 
with perhaps the most common if not the most powerful 
being to guess the solution. (This can be turned into a 
quadratic equation for /i if one expands £>(/i) to lead- 
ing order in L2/L0. There is a small loss of accuracy in 
such an expansion, which is not compensated by the even 
smaller gain of an explicit analytic solution.) This result 
provides a way of measuring the friction coefficient. 



With 7 having been obtained from the measured con- 
tact slope and the calculated rate of change of tip position 
with angle, Eq. (|21)1) . one can now give the surface force 
as a function of separation. From Eq. (fT5|). the angle 
deflection is 



t (z p )= 1 - 1 V cxt (z p ) + E» 



: TTicxt 



(z 



pi- 



rn 



Here the contribution of the linear extrapolation of the 
asymptote of the surface force, F^ xt (z p ), which is as- 
sumed a known function and which was removed from 
the raw voltage signal, has been added back to give the 
full deformation angle. Note that it is the non-contact 
value of the conversion factor, E nc = E(n = 0) that is 
used here. Inserting this into Eq. (|T8|) gives the surface 
force F GXt (z p ). Explicitly in terms of the voltage it is 



F oxt (z p ) 



^V°*\z p ) + F^(z p ), h>0, 

V t t ( 2 9) 

F oxt (zp)+F b cxt (z p ), h = 0. 



S ext 7 

From Eq. (|2T)|) the separation is 



h ext (z p ) = z p + 



Zp + 



(30) 



1 fP nc 



The constant Zg xt is calculated so that h ext = when the 
surfaces first come into contact (see next). The separa- 
tion equation is normally used explicitly for h > 0. These 
three equations are written explicitly for extension; for 
retraction in contact change the superscript 'ext' to 'ret', 
including £ cxt = E{ji) =*> E lct = E(-fi). 



5. Zero of Separation 

The constant Zg xt remains to be determined. The con- 
ventional way of establishing the zero of separation is by 
eye, which is to say the force curve is shifted horizontally 
until it looks 'right'. The problem with this is that there 
is often ambiguities in identifying first contact, particu- 
larly when one has a steeply repulsive surface force prior 
to contact. Also the constant that gives h = at first 
contact may be different to the constant that gives h = 
for most of the contact region, even when the correct 
friction coefficient is used, as will be demonstrated by 
explicit data below. Finally, choosing the zero of separa- 
tion by eye introduces a psychological element into the 
analysis and the potential for personal bias that would 
be best removed by having a mathematical algorithm for 
finding contact. 

It may be objected that identifying the contact and the 
base-line regions already introduce some form of psycho- 
logical bias into the analysis of the experimental data. 
However, it turns out that the various fits are not very 
sensitive to the choice of the region used for the fit, within 
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reason, and the result do not vary significantly with dif- 
ferent choices. The zero of separation, however, feeds 
directly into the final result, and a difference of as small 
as 1 nm can quantitatively effect the values of parame- 
ters that one is trying to measure (e.g. the slip length 
in drainage flow can be of the same order), and it can 
even qualitatively effect the physical interpretation of the 
data. 

The strategy is to define zg xt so that the separation is 
exactly zero at first contact, (and to define zj et so that 
h = at last contact). 'First' (or 'last') contact is defined 
to mean the point at which the extrapolated base-line 
voltage intersects the extrapolated contact voltage. This 
definition is precise and unambiguous, it is able to be 
calculated mathematically, and it is physically reasonable 
and in accord with one intuitive understanding of the 
meaning of contact. 

It should be understood that there are conceptual 
problems with the meaning of 'separation' at the molecu- 
lar level. The separation as defined here, h = z p + z t + zq, 
is not precisely zero over the whole contact region. It 
rather measures the difference between changes in the 
piezo-drive position and changes in the tip position. It is 
positive when there is a protuberance on the substrate, 
and it is negative when there is a depression. It is also 
negative when compression of a deformable surface oc- 
curs. Hence the separation h in contact really gives a 
topographic map of the substrate. The zero plane of 
the map is here defined as the plane passing through the 
point of first or last contact. 

Let Zp Xt be the piezo-drive position at first contact, 
and let V° xt = V cxt (Zp Xt ) be the corrected voltage at 
first contact. (In the linear case, it makes no difference 
to the results what point is selected for z^.) In con- 
tact, V cxt (z p ) = V c cxt + " z pf}- The position at 
which the voltage in contact extrapolates to zero, which 
is defined as first contact, is 



ext 
''pcb 



K ext /U ext . 



(31) 



(The base-line corrected voltage is zero, and so this is the 
same as the intersection of the base-line and contact raw 
voltages.) When the voltage is zero the angular deflection 
is 



next / ^ext 



/ ^ext \ pne poxt / ext \ 

\ Z pcb) — & r b \ Z pchl- 



(32) 



Inserting this into the equation for the separation, and 
setting the latter to zero, ^ cxt ( z p X t,) = 0' gives the shift 
constant, 



-'peb 



iriext / _ext \ 
~ r b \ Z pcb)- 



(33) 



6. Effective Spring Constant 

The conventional modeling of the atomic force micro- 
scope is not only linear but also effectively takes the can- 
tilever to be horizontal. Ignoring the tilt is equivalent to 



equating the cantilever deflection to the vertical position 
of the tip, x = z t . The relationship between the mea- 
sured photo-diode voltage and the vertical tip position is 
given by the calibration factor obtained from the slope 
of the contact region. In this case, the effective spring 
constant that gives the non-contact force is 



fc cff = -i = a nc /£ B 

Zt 



(34) 



The constants E nc = E(ii = 0) and a Ilc = a{n = 0) are 
defined in Eqs (|18l) and (l20|) . respectively. 

The difference between the cantilever spring constant 
ko = 3B/Lq and the effective spring constant k e s can 
be substantial. For the case analyzed in detail below 
(Lq = 110 /zm, R = 10.1 /im, 9 = —11°), the cantilever 
spring constant is fco = 1.37 N/m and the effective spring 
constant is k e g = 1.68 N/m. 

In using the equations for the tilted cantilever to con- 
vert measured atomic force microscope data to force, 
one should use the cantilever spring constant. In using 
the equations for the horizontal cantilever (simple spring 
model, the conventional approach) to convert measured 
atomic force microscope data to force, one should use 
the effective spring constant. In calculating a theoreti- 
cal force curve modeled with the cantilever as a simple 
spring, one should also use the effective spring constant. 

Finally, in Ref . [13 the correct equations for the thermal 
calibration of the atomic force microscope cantilever were 
given. In that paper the cantilever spring constant was 
denoted ko (here also denoted fco), and the effective force 
measuring spring constant was denoted k (here denoted 
fc e ff) and was given in terms of the cantilever spring con- 
stant in Eq. (17) of Ref. [T^- In the present notation, the 
correct thermal calibration method gives the cantilever 
spring constant as 



fco 



— 2/3 c b^O 



3[C + 2L 2 S /L 
x 0.7830 



(35) 



6fc B T 



[Cl + (3L 2 S C /L ) + 3S 2 L 2 2 /Ll] 



Here the measured quantities are /? c b = AV/Az pi which 
is the contact slope evaluated near the base-line voltage, 
(the average of the extend and retract values), /r, which 
is the resonance frequency of the first mode in Hz, PdCj 
which is the direct current power response in V 2 Hz _1 , 
and Q, which is the quality factor. The length of the rigid 
part at the end of the cantilever, L±, defined in earlier 
analyse o 10 ' 14 i 15 has here and throughout been set to zero. 

This inserted into Eq. (I3"4l gives the effective spring 
constant for use when the cantilever is modeled as a sim- 
ple spring (i.e. tilt neglected) , which is usually the case in 
the linear analysis of measured data and the theoretical 
modeling of force-separation curves. 
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7. Effective Drag Length 

The above procedure for analyzing the measured data 
removes the constant force due to the drag on the can- 
tilever from the extension data, and its equal and op- 
posite value from the retraction data. In some case it 
is useful to have available an explicit value for this drag 
force. 

Like the drainage force, and unlike the virtual deflec- 
tion, the contribution to the gradient of the base-line due 
to thermal drift is equal and opposite on extension and 
retraction. 

With (3 D xt and j3 et being the measured base-line slopes 
of the raw voltage as defined above, and F^ xt = — F£ et 
being the gradient of the drainage force in the base-line 
region, then the gradient of the voltage due to thermal 
drift is 



cxt 
b 



PI* 



TTiCXt 



Eiret 
^b 



(36) 



This is equal and opposite to the gradient on retraction. 
With z Pi t U rn being the turn point at the end of the ex- 
tend branch and the beginning of the retract branch, it is 
readily shown that the constant drag force on extension 



pcxt 
drag 



rpcxt 



1 



2jE n 



2-/E nc dz F 



f/cxt 



[2-V 



yret 



(37) 



5 pb 



r ret] 
Z pbJ 



One can define an effective drag length from 



= -6^i p cxt L cff . 



(38) 



One should not take L e ff too literally, but it is expected 
to be somewhat less than the length of the cantilever, 
typically one third to one half of Lq. It should be in- 
dependent of the drive velocity, although because it is 
derived from the difference in the base-line values, it can 
have relative large errors, on the order of 10-20% (see 
results below). 



a spherical probe the dependence on angle is practically 
negligible, and L 2 (9 to t) can be replaced by L2, or even 
by 2R. 

The vertical position of the tip depends upon the de- 
flection and the deflection angle. Again simple geometry 
shows that 



Zt 



J0tot + ^(tftot) sin(0 to t)0tot - L 2 S e . (40) 



Since it is the change in tip position that is important, 
this has been chosen to be zero in the non-deflected state. 
This equation is the major source of non-linearity in the 
cantilever. 

The force F and torque r on the cantilever, which were 
treated in Ejll Al are a linearly proportional to the normal 
force F z and to the lateral force F y , with the proportion- 
ality constant depending upon the lever arm and the tilt 
angle, 



F = F z cos 



and 



FzLo 



t) sin0 t . 



hot + F v siaOtat, 



FyL 2 (0 t ot) cos#tot- 



(41) 



(42) 



Using the linear model of friction, F y = ±fiF z , ijHB 21 
the force and torque are linearly proportional to the sur- 
face force F z . For extension in contact one has 



F = [cos 6»tot + fi sin tot ] F z 



and 



= L 2 



[sin 9 tot - fi cos 6 tot] F z . 



(43) 



(44) 



(Of course for retraction in contact, fx =£- —fj,, and out 
of contact, \i = 0.) Inserting these into the standard 
cantilever equations (J7| and ([8j gives the deflection and 
deflection angle as linearly proportional to the normal 
force (on extension in contact), 



T 3 
3B 



coat 



L 2 L 2 ( 



2B 



/i sin ( 



suit 



^cos6» tot ] > F z , (45) 



III. NON-LINEAR CANTILEVER OR 
PHOTO-DIODE 

A. Non-Linear Cantilever 

For a spherical colloid probe of radius R, simple geom- 
etry, Fig. [51 gives the lever arm as 



and 



= Ry/2 + 2cos( 



(39) 



be written L 2 = 
For the case of a tipped can- 



In the undeflccted state this will 
L 2 (0 O ) = i?V2 + 2cos0 o . 
tilever, L 2 (6 to t) is equal to the length of the tip, and there 
is no dependence on the deflection angle. In fact, even for 



f L 2 

= < ^ [cos 6 tot + M sin ^tot] 



L0L2 (Otat ) 



B 



[sin e to t - A* cos 6» to t] | F z . (46) 



Because the total angle depends upon the deflection an- 
gle, f?tot = $0+^1 these represent a non-linear relationship 
between the force, deflection, and deflection angle. This 
last equation is best written by taking the proportional- 
ity function over to the other side, which gives the force 
as an explicit non-linear function of the deflection angle, 

F z {e). 
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B. Non-Linear Analysis of Measured Data 



(possibly an extrapolation beyond the region of the fit), 



The non-linear analysis in this section holds for both 
sources of non-linearity: the cantilever treated explicitly 
in the preceding subsection and the photo-diode non- 
linearity, for which no specific model is given. In later 
sub-sections, one or other of these will be turned off. 

The raw measured voltage V is a non-linear function 
of the effective total angle, 8 to t = 6 tot + &b — 0m = @o + 
8 + 8 — 8 D {. Note the distinction between the physical or 
relevant total angle dtot and the effective total angle #tot- 
The physical contributions are the tilt angle 8q, which 
is known, and the deflection angle due to the surface 
force, 8, which is to be obtained. The base-line angle 8^ 
is essentially an artifact arising from thermal drift and 
virtual deflection, and it includes the angle due to drag 
force, here assumed constant (but see Refs. [l2lfl3l where 
variable drag is shown to occur for weak cantilevers), 
and the angle deflection due to the linearly extrapolated 
asymptote of the surface force, #bf- In the non- linear 
case, whether it be the non-linear cantilever or the non- 
linear photo-diode, one cannot just subtract the base-line 
voltage from the raw voltage to obtain a corrected voltage 
that gives 8 directly. However, the total angle is a linear 
function its component parts, and so the immediate tasks 
are to obtain the effective total angle from the measured 
raw voltage, 8 t ot(V), and to obtain the base-line angle as 
a function of the piezo-drive position, 8h(z p ). 

Do a non-linear fit of the measured photo-diode voltage 
in the contact region on extension, 



z c p f(V)=a c * t + aT t V + 



2 C X ty2 



tyi 



(47) 



It is best not to use too many terms in the fit. Also 
one should be certain that the fit begins just after first 
contact, and that it ends before any anomalies associated 
with the turn around point at the end of the extension 
branch. 

In the non-linear cantilever case, one regards the an- 
gular deflection 8 as the independent variable. From the 
non-linear Eq. (j46|) . one can calculate the surface force 
F z (8) , and inserting this into the non- linear Eq. (|45p one 
can calculate the deflection x(8). Inserting these into the 
non-linear equation Eq. ([40)) one can calculate the tip 
position zt{8). Hence from these non-linear cantilever 
equations, Zt{8;{i) is easily calculated. A non-linear fit 
can be made to this, (in contact and on extension, /i > 0), 



9^(zt) 



bTzt 



bT^t +&l xt ~ 3 



zt + 



(48) 



The leading coefficient is b cxt — a cxt , which was given 
analytically in Eq. ([20]). One has 8^l c (z t ) = 8 +8 c c xt (zt). 

As in the linear case, the measured voltage is fitted to 
a straight line in the base-line region, Eq. (f2~Tj) . 



V b ex \z p ) = V b ext + ft 



?pbj 



(49) 



Now the raw contact slope from the non-linear fit 
Eq. (jJTJ), evaluated at the base-line voltage constant, V h ext 



pcxt 
Pch 



AV" 



2a cxtye X t 



3a cxt 



+ 



(50) 



The tip position is a linear function of the deflection 
angle in the base-line region, and the gradient is 



d8_ 

dz t 



(51) 



2 t =0 



D^Co + L 2 S + L 2 C 8 - R 2 S$8 Q /L 2 ' 

where D nc = D(0) = [2L^C + 3LlL 2 S ]/[3LlC + 
6L L 2 S ], as follow from Eqs ^ and (120)) . 

The raw voltage is the same function of the total angle 
in and out of contact, and on extension and on retraction. 
Hence the rate of change of raw voltage with total angle 
(equivalently, deflection angle) in the base-line region, 

base-line voltage. One has 



7r 



AV 
A8 

Az t 
A8 



(52) 



:b,cxt 



Az v . 



:b,ext 



Az t 



AV 



cb.cxt ^ Z P 



cb.cxt 



ocxt l „ ext 



(53) 



Here a^ xt = 



a-hip) is given by Eq. (|20)) with \i > 0. One 
can similarly obtain 7^ ot . The correct value of /i makes 
7t> xt = 7b 0t - Actually, the correct /x should make AV j A8 
equal on extend and retract for all voltages. 

With this conversion factor, the angle corresponding 
to the base-line voltage is 



next 
?b 



(z p ) = E^Fl 



( 7 r)-^b xt fe 



J pb J- 



(54) 



This makes the base-line angle at z^ equal to that due 



to the surface force alone; E DC = E(0), Eq. (fl8|). con- 
verts force to angle in the linear regime of low force. The 
change in angle from its value at zj^ xt , includes the con- 
tribution from the linear extrapolation of the asymptote 
of the surface force. The latter is 



8 c h f(z p ) = E nc {F^ t + F^ t '[ 



- p * 



(55) 



not just 



These last two equations apply for all 
in the base-line region. The constant factors do not 
change value in or out of contact. One can see that 
^b Xt ( z p) — ^b Xt ( z p) i s the virtual angle without any con- 
tribution from the surface force asymptote. Hence #tot = 
8 + - 6»^ xt (z p ) - 8lf(z p ) is the total effective angle 
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(physical plus virtual) without double counting the sur- 
face force. 

From Eq. (|4"?| . z pc (V) is known in contact. The tip 
position is related to the piezo-drive position in contact 
by = h = Zr, 



-pc 



+ z tc + zZ xt , or 



Ztc(V) 



(V). 



(56) 



(Strictly, this holds at one particular position in contact. 
See the discussion in the third and fourth paragraphs of 
< 3II B 51 ) The constant ZQ Xt will be determined explicitly 
in the following subsection; in the mean time it can be 
regarded as an arbitrary horizontal shift that establishes 
contact at zero separation. Using successively Eqs (|47| . 

and (l4"8"l) . one can calculate the deflection angle for 
any given raw voltage in contact, 9 C (V). In contact, the 
total effective angle as a function of raw voltage is 

Ototdv) = e + 9?\v) + e^( Zpc (v)) - e° h f(z pc (v)). 

(57) 

This is the desired relationship between the raw photo- 
diode voltage and the total cantilever angle. It is most 
convenient to do a non-linear fit of this (necessarily for 
voltages in the contact region) 



gcxt 
7 tot,c 



(V) 



V- 



iy-2 



(58) 



Although this is explicitly derived for contact on exten- 
sion, it ought equal the analogous result for contact on 
retraction, (assuming the correct value of fi, and the va- 
lidity of the linear friction law) . Either expression can be 
used unchanged for any V measured out of contact and 
on extension or retraction. Hence the left hand side will 
be written simply 9tot(V). 

At a given measured datum {z p , V}, in or out of con- 
tact, on extension, the deflection angle is 



6(z p ,V) = 6 tot (V) - O 



r^(z p ) + e^f(z p ). 



(59) 



The total physical angle is of course 9 to t{z p ,V) — 9 + 
9(z p ,V). From this the force F(z p ,V) follows from 
Eq. (|4"6")l . and the separation h(z p , V) — z p + z t (z p , V) + 
zjf* follows from Eqs (gUJ) and (jfS)) . 

Note that in the non-linear analysis, this equation for 
the separation is applied both in and out of contact; the 
separation is not simply set to zero in contact (c.f. Fig. [I}. 
Non-zero values of h in the contact region give informa- 
tion about the topography and other physical attributes 
of the surfaces, as will be shown below. 



The constant . 



1. Zero of Separation 



remains to be determined in the non- 



linear case. As discussed in the linear case in cjHB 51 
above, there are good reasons for desiring a mathemati- 
cal algorithm for determining the zero of separation that 
avoids guess-work or personal bias. As in the linear 



case, one defines the plane of zero separation as pass- 
ing through the point of first (or last contact), and one 
defines the point of first (or last contact) as the point of 
intersection of the extrapolations of the contact and of 
the base-line raw voltages. In the non-linear case this is 
accomplished as follows. 

Begin by defining z™^ as the piezo-drive position at 
which the fitted voltage in contact, Eq. (jT7|) . equals the 
constant part of the base-line voltage, V b ext , 



-*pcb 



ext 

a 



(60) 



Now a small correction to this will be made such that 
^pcb* 1S * ne piezo-drive position at which the fitted volt- 
age in contact equals the actual base-line voltage at that 
position, viz. z^*'* = Usim ? a Taylor 

expansion to linear order about V^ xt , one has 



ext,* 
Z pcb 



^pcb 



d ycxt 



T TCXt 



/ ext,*\ 
( Z pcb ) 



or 



Scxt „cxt 
Pcb z pcb 



ext 
1) 



-pcb 



T/CXt 



z pb 



ocxt ~cxt 
Pb z pb 



-pcb 



/OCXt 

Pcb 



OCXt 

Pb 



(61) 



(62) 



Since in general 



A' 



ox t 
cb 



> 



the difference between 



cb 



and z™^ is generally small, in many cases negligi- 



At this particular position, the voltages extrapolated 
from contact and from the base-line are equal, V* = 
^c( z pcb*) = ^b( z pcb*)' an d hence the corresponding ef- 
fective total angles must also be equal, 9^ C (V*) = 
9^1 h (V*). The left hand side is 



z {V*)=9 + 9r{V*)+9r{z c p T) 



next/ ext,*\ 
7 bf \ z pcb ) 



and the right hand side is 



Cb(v* 



(63) 



(64) 



Hence at this particular position, the deflection of the 
cantilever due to surface forces is equal to the linearly 
extrapolated deflection due to the base-line surface force, 



next/ cxt,* \ 
7 bf l^pcb ) 



(65) 



with #£ xt (z p ) being given by Eq. ([53]). In this linear part 
of the curve, the position of the tip at this angle is 



cxt,* 
z tcb 



next/ ext,* n / nc 
'bf l^pcb )l a b 



E nc I Fu 



Z?ext' r cxt,* 
^b i Z pcb 



^pb 



]}AC,(66) 



with a^ c = a(0) being given by Eq. (I5ip . which is just 
Eq. (E0D with [i = 0, 
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-100 -50 
Piezo-drive (|im) 



FIG. 3: Calculated deflection angle versus the piezo-drive 
displacement in contact on extension (Lo = 235 /im, R = 
10.11 fim, Go — —11°, and [i — 0.375). The solid curve is the 
non- linear calculation, and the dashed line is the linear cal- 
culation, which gives the tangent at first contact. Note that 
the total drive distance is rather large . 



This particular point is first contact (on extension; it 
is last contact on retraction), and the separation is set to 
zero, h = 0. This gives the off-set as 



-pcb 



6 tcb 



(67) 



This quantity z% xt is used in Eqs ([5d| to (j59")l . which 
equations were not used in its derivation. Any difference 
between zg xt and Zg et mainly reflects the different sign of 
the friction force on extension and retraction, Zt{V, fx) ^ 
zt{V , — fx). The two values bring the non-zero values of 
the separation in contact on the two branches more or 
less into coincidence, as will be demonstrated with the 
data below. 



C. Non-Linear Cantilever, Linear Photo-diode 

Figure [3] shows the deflection angle versus piezo-drive 
position calculated with the non-linear cantilever equa- 
tions given above. In the displayed data, the photo-diode 
does not enter, and so the figure gives an indication of 
the non-linearity due to the cantilever alone. The lin- 
ear result, the straight line, is just the tangent to the 
curve at first contact. One can see that the initial non- 
linear effect is for the angular deflection to increase at a 
slightly decreased rate (i.e. it dips below the tangent; the 
non-linear curve is initially concave down) . Subsequently 
the angular deflection increases at an increasing rate and 
rises above the tangent, (the non-linear curve becomes 
concave up). 

It should be emphasized that the drive distances in 
contact used in Fig. [3l although experimentally realiz- 
able, are on the order of 100 times larger than are used 
in a typical atomic force microscope force measurement. 
The point of the figure is to show what one might qualita- 



tively expect from cantilever non-linearities in an experi- 
ment with large drive distances and a linear photo-diode. 

The non-linearity may be characterised by the horizon- 
tal displacement of the tangent from the non-linear curve 
(the dashed line and the solid curve in Fig. Q} at a given 
voltage. Since the deduced change in piezo-drive position 
is the same as the change in tip position in contact, any 
such horizontal displacement appears directly as a change 
in separation. For the measured data in Fig.[TJ at the end 
of the 300 nm piezo-drive range in contact, the difference 
between the actual piezo-drive and the piezo-drive po- 
sition at which the linear tangent at contact gives the 
same voltage is 54 nm. For the calculated data for the 
non-linear cantilever used in Fig. [3J but over the same 
range of 300 nm, the difference is 1.2 nm. 

One can conclude from these calculations that the non- 
linearity in the photo-diode is about 50 times greater 
than the non-linearity of the cantilever in this series of 
experiments. 



D. Linear Cantilever, Non-Linear Photo-diode 

The present series of experiments are likely representa- 
tive of the norm, and the dominance of the photo-diode 
non-linearity over the cantilever non-linearity is likely the 
rule rather than the exception. Hence it is worthwhile to 
give explicitly the equations required to analyze atomic 
force microscope force data in the case of a non-linear 
photo-diode and a linear cantilever. These simplify the 
full analysis given earlier in this section. Essentially one 
takes the linear cantilever results from mi B 31 and com- 
bines them with the non-linear data analysis results of 

mrsi 

For the linear cantilever, the lever arm is the con- 
stant L2 = R\/2 + 2 cos 6q . Also, the deflection angle, 
force, and vertical tip position are all proportional to 
each other, Eqs QBH2D). 



D ext (9, and z t = 6/a c 



(68) 



The superscript 'ext' (fj, > 0), is replaced by 'ret' (/x < 0), 
or by 'nc' (p = 0) in the respective regimes. 

For the non-linear analysis of the measured data, the 
fit to the measured raw voltage is as in Eq. (|47|) , 



z;f(V) = a + 



cxt 



V 



'[V] 3 



(69) 



As in the linear case, the measured voltage is fitted to a 
straight line in the base-line region, Eq. (|2"T|) . 



v h ext (z p ) = v b ext 



z P b}. 



(70) 



Now the raw contact slope from the non-linear fit eval- 
uated at the base-line voltage constant, V£ xt is given by 
Eq. pi, 



Scxt 
Pcb 
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3a| xt [F b ext ] 2 



(71) 



The rate of change of angle with tip position in the base- 
line region, the a£ c of Eq. (|5Tj) , is the same as the a(0) 
of Eq. (|20|) . The conversion factor for voltage to angle 
in the base-line region given by Eq. ([52]) is unchanged, 



7 ext 



-/3 c °£V«b xt - 




The base-line angle and the base-line angle due to the 
surface force are unchanged from Eqs (|54)l and (1551) . 
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\z p ) = E nc F£ + ( 7 D 



ext \ — 1 /ocxt r 



Pb 



J pb J; 



(72) 
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FIG. 4: Measured^ and calculated force versus separation for 



#b Xt (z ) = E nc \F£ xt + Fu Xt ' [z — z cx ^]\ . (73) a drive velocity of i p = 2 /im s 1 , viscosity r] = 51.62 mPas 
' " 1 I b P pb J friction coefficient a = 0.35. and cantilever snrine const; 



The total effective angle in contact is formally the same 
as Eq. JSSJ), 



W = *° + ^ Xt (^) + ^ Xt (^ P c(^)) - flbf^peOO), 

. (74) 

with the deflection angle now being linear in the tip po- 
sition, 



friction coefficient fi = 0.35, and cantilever spring constant 
fco = 1.37 N/m. The symbols are the atomic force micro- 
scope data (upper is for extension, lower is for retraction), 
analyzed with the non-linear photo-diode, linear cantilever 
algorithm. The almost overlapping dashed curves are the cal- 
culated drainage force with a slip length of b = 3nm and 
effective spring constant of k e g = 1.68 N/m. The right inset 
magnifies the measured data in contact (jagged solid curve is 
extension, jagged dashed curve is retraction). 



t {V)=a c 



(75) 



At a given measured datum {z p , V}, in or out of contact, 
on extension, the deflection angle is formally the same as 
Eq. (US]), 



0(z p , V) = 9 tot (V) -6 - 9l xt (z p ) + 0&(zp). (76) 

From this, the linear equations ([55)1 give the force and 
also the separation, h = z p + z t + z$ xt . 

The equations that give the constant that fixes the 
zero of separation are essentially unchanged from mil B II 
One has the crude estimate of the piezo-drive position 
where the base-line intercepts the contact curve, 



cxt ~cxt I cxtT/cxt | cxt [T/ext"|2 I ext [T/extl ci , 

Z pcb — a + a l V b + a 2 l V b \ + a 3 l V b I + 



and the refined version, Eq. (|6"2")l . 



(77) 
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FIG. 5: Magnification of previous figure at small separations. 
The symbols and dashed curves are as in the preceding figure, 
and the dotted curves are the calculated stick drainage force 
(6 = 0nm, k cS = 1.68 N/m). 
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Equation flM) again follows, 6 cxt {V*) = 0gf (z^*), 
which in the present case can be rearranged as 



r ext _ ext,* 
Z — Z pcb ' 



[E nc /a cxt ] {Fg* + Fg*'[ 



ext' r_ ext >* _ 

z pcb ''pb 



£h"J } ■ 

(79) 

This exactly the same as given in mil B 11 as one might 
expect since linear analysis suffices for the base-line. 



IV. RESULTS 

Figure @] shows the force versus separation data that 
corresponds to the raw measured data presented in Fig. [I] 
These and all of the experimental results presented in this 
section are analysed using the non-linear photo-diode, 
linear cantilever algorithm presented in the preceding 
subsection, unless explicitly stated otherwise. Four co- 
efficients were used in each of the non-linear fits, do, oti, 
0,2, and 0.3 on each branch. The measurements in Fig. 0] 
were performed at low velocity, z p — 2 /ims -1 , and so the 
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magnitude of the drainage force is comparatively weak. 

The figure includes the calculated drainage force, 
which is given by the Taylor equation with slip, 



-6m]R 2 h 



f(h), 



where the slip factor isi£ 
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with 6 being the slip length. The case 6 = (equivalently, 
f(h) = 1) corresponds to stick boundary conditions. 

A weak van der Waals force (Hamaker constant 1.3 x 
10~ 21 J m~ 2 ) and short range repulsion (equilibrium sep- 
aration 0.53 nm) derived from a Lennard- Jones 6-12 po- 
tential was included in order to calculate results in con- 
tact. These have no effect for separations greater than 
about 1 nm. Details of the algorithm used to calculate 
the drainage force may be found in Ref. |12| . 

Comprehensive atomic force microscope measurements 
of the drainage force and of the slip length have been 
given in earlier work by the present authors ^ 12 ' 13 albeit 
with the linear analysis. The focus in the present work is 
not directly on the drainage force or the slip length, but 
rather on those aspects of the measurements where the 
non-linear analysis makes a difference. 

Accordingly, one can briefly observe that the agree- 
ment between theory and measurement in Fig. [4] indi- 
cates the overall validity of the non-linear procedures 
developed here. The cantilever spring constant, fco = 
1.37 N/m, was obtained by fitting by eye the data and 
the calculated force at larger separations where the stick 
and slip theories overlap. In fact, at this velocity, 1.3 < 
fco < 1.6 N/m gave equally good fits. At the higher ve- 
locity z p = 20/ims~ 1 , 1.3 < fco < 1.5 N/m were ac- 
ceptable, and at the still higher velocity z p = 50 /ims -1 , 
1.33 < fc < 1.42 N/m fitted the data. (The higher 
velocities gave larger drainage forces and the fits were 
more definitive. These cases are discussed in detail be- 
low.) Based on these fits, the value of the cantilever 
spring constant is taken to be fco = 1.37 N/m, here 
and below. It is estimated that the value has been ob- 
tained with a precision of about ±0.1 N/m. The effective 
spring constant corresponding to this is given by Eq. (f3~4"|) , 
fc c ff = 1.68 N/m. This effective spring constant should be 
used in the conventional linear analysis of the data, and 
in the theoretical calculations that use a simple spring 
model. 

Figure [5] shows the quality of the fit used to obtain the 
cantilever spring constant at the lowest drive velocity. It 
can be seen that a value of fco = 1.37 N/m (equivalently, 
fc c ff = 1.68 N/m) slightly overestimates the magnitude of 
the deflection, more noticeably on the retraction branch, 
and that a slightly higher value would give a better fit 
in this case. Note that in order to reduce the fitting pa- 
rameters from two to one, only the stick theory is used 
for the fit. (In any case the slip theory is virtually coin- 
cident at these large separations.) Also included is the 
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FIG. 6: Fit for the spring constant, (data corresponding to 
Fig- SI- The vertical tip position, zt = F/k e &, is shown. The 
plus symbols are the non-linear analysis of the raw experimen- 
tal data,— ^ the cross symbols are the original linear analysis 
of the same data,— the solid curve is the position calculated 
from the stick drainage force (b = Onm, fe e ff = 1.68 N/m), 
and the dashed curve is the position calculated from the stick 
drainage force (b = Onm, k e g = 1.5 N/m). 



original linear analysis of the same raw data^ and the cal- 
culated stick drainage force that was used to fit a value 
of fe e ff = 1.5 N/m (equivalently, feo = 1.23 N/m). First 
it may be observed that there is only a minor difference 
between the linearly and the non-linearly analyzed exper- 
imental data; the linear data slightly overestimates the 
magnitude of the deflection, and this effect increases as 
the magnitude of the force increases. Second, the value of 
the spring constant used originally is about 10% smaller 
than that fitted here. This gives a slightly better fit to 
the linearly analyzed extend data, but a worse fit on the 
retract curve, which was not taken into account in Ref. 0. 
At smaller separations than shown, the original and the 
present fits slightly improve, but in this region the slip 
length begins to have an effect and so it is better not 
to include it in the determination of the spring constant. 
Although the discrepancy in the spring constant is small 
(errors in other methods of cantilever calibration typi- 
cally exceed 20%), it does lead to a large discrepancy in 
the fitted slip length and to qualitatively different be- 
havior of the slip length with shear rate to that found 
here. 

The overlap between the measured data and the slip 
calculations at short range in Figs @] and [3] indicates that 
the correct slip length, 6 = 3 nm, has been obtained. The 
agreement is surprisingly good for the retract data, which 
was not analyzed in Refs [6l ll2lll3[ but which turns out to 
be very sensitive to the slip length. For example, the 
adhesive force, which is defined to be the minimum of 
the retract force, in Fig. [5] is measured to be — 16.6nN, 
and it is calculated to be —19.1 nN for 6 = 3nm For the 
case of stick (6 = Onm), it is calculated to be — 36.3nN, 
which significantly overestimates the adhesion. The ex- 
tend data is less sensitive to the slip length. For example, 



14 



at h = 10 nm, the measured force is 10.5 nN, the calcu- 
lated slip force for b = 3nm is 10.1 nN, and the calcu- 
lated stick force, b = Onm, is 13.2 nN. Based on this and 
the higher velocity results shown below, it is estimated 
that the slip length has been obtained with a precision 
of about ±0.5 nm. 

The linear analysis (not shown) of the measured ex- 
tend data using the slope of the tangent at first contact 
gives 10.8 nN at h = 10 nm, and, using the slope of the 
tangent at final contact it gives — 16.5nN for the adhe- 
sion. Hence the error due to the linear analysis is quite 
small at these points at this lowest velocity. It should 
be noted that the linear results were obtained using the 
effective spring constant k e s — 1.68 N/m rather than the 
cantilever spring constant, ko = 1.37 N/m. 

It is worth mentioning that the drainage measurements 
reported in Ref. H found a low shear rate limiting slip 
length of bo = 10 nm for the Si-Si case. This is signif- 
icantly larger than is found here. In that investigation 
it was found that the slip theory progressively under- 
estimated the measured force at small separations and 
overall it was a noticeably worse fit than found here. As 
mentioned above, the retract data was not included in 
the fit in Ref. H|. The most likely reason for the larger 
slip length and the degradation in the fit at small sepa- 
rations in Ref. |6| compared to here is the difference in the 
value of the spring constant used: here k e g = 1.68 N/m, 
compared to k c g — 1.5 N/m in Ref.S A detailed discus- 
sion of how the linear analysis led to the contradiction 
between these earlier results and the present ones is given 
in the concluding section. 

The friction coefficient was determined from the mea- 
sured data in Fig. [4] by equalizing the extension and 
retraction calibration factors using the Guided Unbi- 
ased Estimate for a Single Solution algorithm, which is 
a rapid and well-used procedure for solving non-linear 
equations. A value of fx = 0.38 gave 7°£* = 6357V/rad 
and 7^* = 6359V/rad. This is a little larger than the 
value determined by Stiernstedt et al.^ who obtained 
fj, = 0.32-0.35 in four independent friction measurements 
(two probes, axial and lateral methods) for a silica probe 
on a silica substrate. In view of this a value of jU = 0.35 
was used to analyse the data in Figs @] and [5] and all of 
the following figures. 

It was found that the value of the friction coeffi- 
cient varied somewhat with the voltage at which 7 was 
evaluated. Obtaining /i by minimising [7 cxt (V ; /i) — 
7 ret (V;/z)] 2 over the contact region also gave a differ- 
ent value. In some case, particularly at higher drive ve- 
locities, unphysical values of fj, were obtained. It was 
concluded that for these particular measurements in this 
high viscosity liquid the non-linear method could not be 
used to obtain the friction coefficient reliably. As is dis- 
cussed further below, the most likely reason for the failure 
of this method of measuring the friction coefficient in the 
present case is the neglect of the variations in the drag 
force, which are most pronounced in contact. 

Interestingly enough, the choice of the friction coeffi- 
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FIG. 7: Measured^ separation versus displacement on sub- 
strate in contact corresponding to Fig. [4] The solid curve is 
extension and the dashed curve is retraction. 

cient (and the choice of the fitting region) had almost 
no effect on the values of the measured non-contact 
forces. This is discussed further below, but in essence, 
the method of analyzing the non-contact data cancels 
the friction contribution whatever it may be. 

The inset of Fig. U shows the force versus separation 
just prior to and in contact. Overall in contact the non- 
linear analysis gives a quite vertical curve that fluctuates 
about h = with apparent noise on the order of ±0.5 nm. 
The verticality of the contact region of the data analysed 
with the non-linear algorithm is much better than that 
obtained with the linear analysis. As can be seen in the 
inset to Fig. [T] the linear analysis gives a systematic error 
in contact of 20-50 nm. 

This noise in the non-linear data is not entirely random 
or unphysical (hence the word 'apparent') because any 
roughness of the surface gives rise to physical fluctuations 
in the apparent separation, as is discussed quantitatively 
shortly. In fact the presentation of the data in Fig. 0] 
exaggerates the apparent noise in the non-linear analysis. 
This is because the friction force is reversed on the retract 
branch compared to the extend branch, and so a given 
surface force F corresponds to different contact positions 
on the substrate. Hence the disagreement between the 
extend and retract traces when presented in the form of 
force versus separation is entirely as one would expect. 

One needs to compare the two branches at the same 
physical position on the substrate. The basis of the anal- 
ysis of the tilted cantilever is that in contact the probe 
slides along the substrate in the axial direction of the 
cantilever, which is the origin of the friction force. The 
displacement from first contact of the probe on the sub- 
strate is simply y — L 2 0. (Actually, the full expression is 
y = S a x - C L 2 6 = [S D cxt / rct - C L 2 }8, but this is not 
used here.) This is used to plot the apparent separation, 
h = z p + z t + zq, for extend and retract in Fig. [7J As 
in the preceding figure, the linear cantilever, non-linear 
photo-diode analysis was used to obtain Zt(9(V , z p )). 

As mentioned in the text (c.f. the fourth paragraph of 
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£|II B 51) , the separation in contact measures the difference 
between changes in the piezo-drive position and changes 
in the tip position. It is positive when there is a protu- 
berance on the substrate, and it is negative when there 
is a depression. That this is really a topographic map of 
the substrate can be seen by the high degree of correla- 
tion between the extend and the retract traces in Fig. [7] 
There are slight depressions, about 0.2 nm deep, at about 
y = 20 nm and y = 38 nm, and there is a protuberance, 
about 0.5 nm high, at about y — 30 nm. Where there 
is no correlation between the extend and retract traces, 
one can attribute the departure from zero to noise that is 
most likely connected with stick-slip motion of the probe 
on the substrate. One can conclude that the substrate is 
really quite smooth and, despite this noise, the method 
does give reliable topographic information with a resolu- 
tion on the order of 0.1 nm. 

It was found that z% xt = —1506. 9 nm and that zj, ct = 
— 1507. 8 nm. These differ by 1 nm, probably due to the 
reversal of the friction force on the two branches, which 
causes a displacement of the vertical position of the tip. 
This 1 nm difference appears to be essential to get the 
very precise overlap on the topographic plot. It is em- 
phasized that the two values of zq were generated by the 
algorithm and no human judgement or adjustment was 
involved. 

The rather pronounced feature in the start of the re- 
tract curve at about y = 42 nm in Fig. [7] is not a phys- 
ical feature. It is an artefact of the reversal of direction 
of the piezo-drive on the change from extend to retract. 
This reversal begins at the start of the retract branch in 
this particular model of the atomic force microscope. At 
this turning point, the probe is essentially stuck at one 
contact position as the friction force reverses direction. 
During this turning phase the present analysis, which as- 
sumes that F y — —/j,F z , is not valid. For the same reason 
the negative separation of about h w — 0.5nm just after 
first contact on extension is likely an artifact of the linear 
friction model. 

Figure [8] shows the measured and calculated drainage 
force for a drive velocity of z p = 20 (ims -1 , which is ten 
times the velocity of the preceding case. The viscosity 
1] = 52.25 mPas -1 has slightly increased due to a tem- 
perature rise in the measurement cell. The cantilever 
spring constant was unchanged, fco = 1.37 N/m. The cal- 
culated slip force (obscured dashed curves) was obtained 
with a simple spring model and the unchanged effective 
spring constant of fc c ff = 1.68 N/m and unchanged slip 
length 6 = 3 nm. The fact that an unchanged spring con- 
stant and slip length fit the measured data equally well 
at this high velocity as at the preceding low velocity in- 
dicates the reality of the values used, the validity of the 
mathematical form for the drainage force, and the reli- 
ability of the non-linear algorithm for the conversion of 
the raw measured signal into quantitative data. 

The fact that the slip length is the same here as for 
the low velocity case (and for the higher velocity of 
i p = 50^ms _1 presented below) suggests that the slip 
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FIG. 8: Measured^ and calculated force versus separa- 
tion for a drive velocity of z p = 20/j,ms _1 and viscosity 
r\ = 52.25 mPas . The symbols are the atomic force mi- 
croscope data (upper is for extension, lower is for retraction), 
analyzed with the non-linear photo-diode, linear cantilever al- 
gorithm with cantilever spring constant ko = 1.37 N/m and 
with friction coefficient of y, = 0.35. The almost completely 
obscured dashed curves are the calculated drainage force with 
a slip length of b — 3nm and k e g = 1.68 N/m. The dotted 
curves are the calculated drainage force for stick boundary 
conditions (6 = Onm and fc e ff = 1.68 N/m). 



length is independent of the shear rate. This conclu- 
sion is reinforced by the fact that the slip theory with 
constant slip length fits the measured data over all sep- 
arations (c.f. Figs [5] and [SJ , since the maximum shear 
rate of the drainage flow increases with decreasing sepa- 
ration. This contrasts with earlier work in which it was 
found that the slip theory increasingly underestimated 
the measured drainage force on extension as the sepa- 
ration approached zeroJ^ As was discussed above, such 
behavior occurs when the spring constant is underesti- 
mated, which was the case in Ref. from which series of 
measurements the raw data used here was taken. It is not 
known at this time whether or not a similar underesti- 
mate occurred in Ref. [l9|. Accordingly, the claim that the 
slip length decreases with increasing shear ratei^ should 
be treated with caution pending a more refined analysis 
of that data. 

The adhesion measured using the non-linear analysis 
is — 82.6nN, and that calculated with the slip length b = 
3nm is — 82.6nN. By way of comparison, the adhesion 
measured using the linear analysis and fc e ff = 1.68 N/m is 
— 80.8nN (not shown), and that calculated with the stick 
theory b = Onm is — 136.6 nN. Note that it is essential 
to use the effective spring constant in order to get such 
reliable linear results. 

The apparent separation in contact as a function of 
position on the substrate is shown in Fig. [§J There is 
a degree of correlation between the extend and retract 
traces, with a noticeable crater appearing in both at 
about y — 20 nm. The amount of noise is noticeably 
larger at this velocity of z p — 20/ims -1 compared to 
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FIG. 9: Measured^' 7 separation versus displacement on sub- 
strate in contact, (corresponding to preceding figure). The 
solid curve is extension and the dashed curve is retraction. 

the data in Fig. [71 which was obtained at a velocity of 
i p = 2 /ims -1 . The physical origin of the increased noise 
is likely that more power is dissipated at the greater ve- 
locity. Two signatures of this change can be noticed. 
First is the reduced correlation between the extend and 
the retract traces. And second is that the amplitude of 
the fluctuations appear larger. The spatial wave length of 
the fluctuations is noticeably larger in Fig. |H] compared to 
Fig. [Jj This indicates that the fluctuations are dynamic 
in nature, and that they probably have the same tem- 
poral frequency in both cases. This point explains why 
the low velocity data is more satisfactory in producing 
a topographic map of the surface: at low velocities the 
high frequency vibrations of the cantilever are averaged 
out in the time it takes to traverse a given spatial feature, 
whereas at high velocities the period of vibration and the 
spatial period are comparable and so they interfere with 
each other. One can conclude that the low velocity data 
gives more reliable results for the contact region than the 
high velocity case. 

This last point is important. The friction coefficient 
used in the non-linear analysis in the high velocity case 
was y, = 0.35, which was taken from the low velocity data. 
(A non-linear fit in the contact region was still performed 
and used for the high velocity case.) This value of // in 
the present case of i p = 20/xms -1 gives calibration co- 
efficients = 5660V/rad and 7™* = 6425V/rad. (In 
the analysis of the data, these values of 7 were used on 
each respective branch.) An unrealistically high value 
of fi = 1.3 is required to bring these into agreement at 
7cb = 6250V/rad. Although the required friction coef- 
ficient is unrealistically high, using it only changed the 
calibration coefficient by about 3%; the measured value 
of the adhesion remained unchanged at — 82.6nN. 

Besides the increased noise due to cantilever vibration, 
there are two further reasons why the high velocity data 
is unreliable in contact. First, in addition to the friction 
force, an axial drag force acts on the probe parallel to the 
substrate as it slides on the substrate. This axial drag 



force is not taken into account in the present analysis. It 
is exacerbated in high viscosity liquids and at high ve- 
locities. Second, the drag force normal to the cantilever, 
which is distinct from the drainage force, is not in fact 
constant as assumed here, but varies with loadJ^H (For 
determining the slip length, it is not necessary to use 
here a variable drag model because the cantilever is rel- 
atively stiff.) 12 i 13 This variable drag effect increases with 
increasing viscosity, increasing velocity, and increasing 
force. Although for stiff cantilevers such as the present 
the variation in drag is relatively negligible in the non- 
contact region, in contact the variation is independent of 
the stiffness of the cantilever and so it can be expected 
to affect the present results. A simple calculation shows 
that here the axial drag force is several orders of mag- 
nitude smaller than the friction force. Hence it is most 
likely that it is the present neglect of the variable drag 
force in contact that makes the determination of the fric- 
tion coefficient in the present series of measurements un- 
reliable. 

The question naturally arises that if the friction co- 
efficient is unreliable for a given velocity, why can one 
use the non-linear contact fits at that velocity? There 
are two answers to this, one pragmatic and one rea- 
soned. First, the evidence is that the fits are reliable 
and produce quantitative measured results, namely that 
the measured adhesion is insensitive to the value of the 
friction coefficient used in the analysis, and the measured 
force agrees quantitatively with the calculated force. Sec- 
ond, the non-linearity in the raw voltage measured in 
contact arises almost entirely from the photo-diode; the 
cantilever, including the friction, contributes to the lin- 
ear part. The two drag forces just mentioned and ne- 
glected here also contribute to the linear part. Similar 
to the friction force, they reverse sign between extend 
and retract. The method of analyzing the non-contact 
data essentially cancels the friction contribution, and so 
it would also have canceled these drag contributions if 
they had been taken into account. This is the justifi- 
cation for using the friction coefficient measured at low 
velocities in combination with the non-linear fits applied 
at the current velocity. 

Figure [TU] shows the measured and calculated force 
for a drive velocity of z p = 50/ims -1 . The measured 
data are analyzed with the non-linear algorithm of mil Dl 
(symbols), and with the linear algorithm of mi Bl us- 
ing the tangent at first (extension) or last (retraction) 
contact (dash-dotted curves). The cantilever spring con- 
stant, ko = 1.37 N/m, slip length, b = 3nm, and friction 
coefficient, /i = 0.35, have been fixed at the values deter- 
mined in the low velocity case, z p = 2 /mis -1 . 

With this value of /.i, in the present case of 
i p = 50/ims -1 the calibration coefficients are 7^' = 
6263V/rad and 7^ = 6416 V/rad. A value of /* = 0.57 
gives 7*jf = 6361 V/rad and Y c b = 6365 V/rad. 

The adhesion measured using the non-linear analysis 
is — 153 nN, that measured using the linear analysis is 
— 147 nN , that calculated with the slip length b = 3 nm 
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FIG. 10: Measured®'- and calculated force versus separa- 
tion for a drive velocity of z p — 50/zms -1 and viscosity 
j) = 52.25 mPas" 1 . The symbols are the atomic force mi- 
croscope data (upper is for extension, lower is for retraction), 
analyzed with the non-linear photo-diode, linear cantilever al- 
gorithm with cantilever spring constant fco = 1.37 N/m and 
with friction coefficient of /i = 0.35. The mainly obscured 
dashed curves are the calculated drainage force with a slip 
length of b — 3nm and k e g — 1.68 N/m. The dotted curves 
are the calculated drainage force for stick boundary conditions 
(b — Onm and fc e ff = 1.68 N/m). The obscured dash-dotted 
curves are measured data analyzed with the conventional lin- 
ear approach using fe e fl = 1.68 N/m, with the zero of separa- 
tion established as described in B 51 

is — 147 nN, and that calculated with the stick theory 
b = Onm is -230 nN. 

In the case of the linear analysis, the zero of separa- 
tion was established as described in All B 51 This removes 
any human intervention or bias, which can be a problem 
in the conventional linear analysis where each curve is 
shifted horizontally so that it gives h = at what ap- 
pears by eye to be first or last contact. In fact, in the 
absence of the non-linear result it would have been diffi- 
cult to establish contact with any confidence in this case 
due to the smooth and continuous nature of the forces 
and the lack of verticality in contact in the linearly an- 
alyzed data. These results with negative gradient and 
negative separation for the linear analysis in contact is a 
significant failing of the approach. 

The verticality of the linear analysis in contact for 
forces F z <E [— 150 nN,150 nN], suggests that the photo- 
diode is linear in this range. Since the non-contact forces 
also lie in this range, one can understand why the linear 
analysis with effective spring constant works in this high 
velocity case. 

It can be noted in Fig. [TU] that the maximum magni- 
tude of the drainage forces out of contact are about on 
the limit of photo-diode linearity. At higher velocities 
they will enter the non-linear regime, and so one would 
expect the linear analysis to become progressively less 
reliable as the velocity is increased. 

It can also be seen that on the extension branch, in this 
figure and in the preceding figures, there is no overlap be- 



FIG. 11: Measured^ separation versus displacement on sub- 
strate in contact corresponding to preceding figure. The solid 
curve is extension and the dashed curve is retraction. 

tween the contact voltages and the non-contact voltages. 
However, the retraction branch in contact encompasses 
all of the non-contact voltages. Whereas extrapolation of 
the contact fit is required on extension, only interpolation 
is required on retraction branch. 

There is a noticeable discontinuity in the non-linear 
measured data at final contact on retraction, and at first 
contact on extension (obscured). This is an artefact of 
the linear friction model, which is zero out of contact, 
and jumps discontinuously to a non-zero value in con- 
tact. Whenever the load is non-zero at such a point, 
the friction force is discontinuous and it gives rise to a 
discontinuity in the analyzed data. 

Figure lTTl shows the apparent separation versus contact 
position on the substrate at this highest velocity. In this 
case it is difficult to see a correlation between the extend 
and the retract traces. The spatial wave length of the 
noise is large, as one might expect of temporal vibrations 
at this high speed, and they appear to dominate any to- 
pographic features. The pronounced feature at the start 
of the retract trace at y = 45 nm is an artefact of the 
turn point and the inapplicability of the linear friction 
model there. Likewise, the pronounced feature near final 
contact, y < Onm on retraction, which can also be seen 
in Fig. [9j is also likely an artifact of the linear friction 
model, F y = —\iF z . In this region on retraction just prior 
to pull-off there is a pronounced adhesion and F z < 0. 
The linear friction model is only valid, if at all, for pos- 
itive loads. In this context it is worth mentioning that 
there is a modified form of Amontons law, F y = fj,[F z —A], 
where A is the adhesion^ There is some experimental 
support for this ; 16 ' 17 but it has not been explored in the 
present investigation. 

In mi B 71 a method of extracting the drag force was 
given, which was characterized by an effective length that 
was expected to be somewhat less than the length of 
the cantilever. The validity of the procedure could be 
checked from the extent to which the effective length was 
independent of the drive velocity. For velocities z = 2, 20, 
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50, and 80 /inns -1 , it was found that L e g = 83.2, 92.9, 
77.7, and 88.1 //m, respectively. These are remarkably 
consistent. There was some dependence on the choice of 
the base-line region. In the worst case, z = 2/ims _1 , 
reducing the range of the base-line by a factor of almost 
3, from 1.83 /im to 0.65 /im increased the drag length by 
about 20%. 



V. CONCLUSION 

This paper has investigated the causes of non-linearity 
in the contact region of atomic force microscope measure- 
ments of surface forces, and has developed an algorithm 
for analyzing data when such curved compliance occurs. 
Two sources of non-linearity — large cantilever deflection 
and photo-diode response — were identified and investi- 
gated. It was concluded that in a typical case the non- 
linear photo-diode response was the dominant contribu- 
tion to the observed non-linearity in contact. 

A relatively simple algorithm for analyzing raw experi- 
mental data was developed. The numerical algorithm in- 
voked a non-linear polynomial fit to the measured voltage 
in the contact region, and was found easy to implement 
with a spread sheet. 

The first advantage of the algorithm is that it elim- 
inates the ambiguity in the choice of the contact slope 
(constant compliance factor) that occurs when one has 
a curved contact region. Even when, or, more precisely, 
especially when, the non-contact forces are in the linear 
response regime, one still needs this calibration factor to 
convert the measured photo-diode voltage to force and 
separation. It significantly improves the quantitative re- 
liability of the atomic force microscope to have an algo- 
rithm that gives it correctly and unambiguously. 

A second advantage is that the algorithm eliminates 
non-physical behavior of the analysed data in the contact 
region that is an artifact of the linear analysis. The re- 
sults show that the non-linear analysis yields reliable to- 
pographic data in the contact region with sub-nanometer 
accuracy. 

The complete non-linear analysis of cantilever deflec- 
tion is useful even though it turns out that in this in- 
stance it is not required in full. The full analysis was 
simplified to the linear case but it still included the ef- 
fects of cantilever tilt, friction in contact, and torque due 
to the extended probe. These are often neglected in the 
conventional analysis but they are required for accurate 
and reliable analysis of raw experimental data. In par- 
ticular, a third feature of the analysis is that it gives 
the relationship between the intrinsic cantilever spring 
constant, which is a material property of the cantilever 
and which is the quantity usually measured by calibra- 
tion procedures, and the effective spring constant, which 
is the quantity required to convert the measured vertical 
cantilever deflection to a surface force. The difference 
between these two can typically be 10% or more, and us- 
ing the wrong one can lead to unacceptable errors in the 



measured surface force. 

A fourth feature of the data analysis algorithm is the 
accounting of measurement artifacts such as virtual de- 
flection, thermal drift, cantilever drag, and long-range 
surface force asymptotes. The proper treatment and 
elimination of these is automated using relatively straight 
forward linear corrections. 

A fifth innovation is an automated numerical proce- 
dure for fixing the zero of separation. The advantage of 
this is that it eliminates the ambiguity that exists in the 
conventional analysis and it avoids human intervention 
and bias. An unambiguous and precise definition of zero 
can be essential when finer details of surface forces in the 
non-contact region are required. 



1. Reconciliation with Earlier Results 

The non-linear algorithm was applied to measured raw 
atomic force microscope data that had previously been 
analyzed using the conventional linear approach. The 
validity of non-linearly analyzed data was confirmed by 
the quantitative agreement with the calculated drainage 
force over a range of velocities^ In particular, quantita- 
tive agreement occurred for the measured drainage adhe- 
sion, which had not previously been analyzed in detail. 

Some effort was made to test the conventional linear 
analysis when the data showed a non-linear contact re- 
gion (curved compliance). It was found that provided 
the correct effective spring constant was used in the lin- 
ear analysis, the measured pre-contact forces agreed with 
the non-linearly analyzed data. Obviously, if the intrinsic 
cantilever spring constant was used instead, as is often 
the case in conventional analysis, the non-contact force 
was underestimated by about 10%. Also, in order to get 
agreement, the zero of separation had to be established 
using the (linear) algorithm given here. 

The slip length obtained on the basis of the present 
non-linear analysis was 3 nm and it was found to be in- 
dependent of shear rate. In the previous work where 
the data was analyzed with the linear algorithm^ the 
slip length was found to decrease with increasing shear 
rate, and to have a low shear rate limiting value of 10 nm. 
These qualitative and quantitative discrepancies were at- 
tributed directly to an underestimate of the cantilever 
spring constant in the previous work^ which is a conse- 
quence of the linear analysis used there. 

In the light of the present results and analysis one 
can explain the reasons for the two major discrepancies 
between the present results and the earlier results: the 
much larger slip length, and the decrease in slip length 
with increasing shear rate found earlier. The earlier anal- 
ysis used as the calibration factor (the conversion fac- 
tor between photodiode voltage and cantilever deflection) 
the gradient of the contact extension curve measured 
at first contact. This corresponds to the lowest exten- 
sion contact force but the greatest extension non-contact 
drainage force. Because the contact curve is concave 
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down due to the non-linearity in the photodiode response 
(see Fig. [T]), this value of the calibration factor under- 
estimates the actual response in the non-contact region. 
Hence, the deflection of the cantilever (voltage divided by 
calibration factor) was overestimated in the earlier work. 
Hence the effective spring constant used in the theory to 
fit this measured deflection at large separations was un- 
derestimated. (This explains the earlier k e g = 1.5N/mj£ 
compared to the present k e s = 1.68 N/m.) However, us- 
ing too soft a cantilever in the theoretical calculations 
of the drainage force overestimates the curvature (rate 
of increase with separation) in the calculated deflection 
compared to the actual curvature at intermediate sep- 
arations. For the given spring constant, the only way 
to reduce the theoretical force so that it agrees with 
the calculated force at intermediate separations is too 
increase the slip length, which causes it to be overesti- 
mated. (This explains the earlier 60 = 10 nni|S compared 
to the present b = 3nm.) However, at small separations, 
where the pre-contact voltage is close to the first contact 
voltage, the calibration factor that was used earlier is 
relatively accurate, and the artificial compensation fac- 
tors (too small a spring constant, too large a slip length) 
that worked at large and intermediate separations are 
not required at small separations. Hence for the given 
spring constant, the overestimated slip length that fit- 
ted at intermediate separations needs to be reduced at 
small separations to fit theory to the analysed measured 
data. Since the shear rate increases with decreasing sep- 
aration, this explains the discrepancy between the earlier 
conclusion that the slip length decreases with increasing 
shear rate^ and the present conclusion that it is con- 



stant. This chain of reasoning likely applies as well to 
other measurements! 12 i 13 ' 19 

Besides substantially improving the reliability of the 
non-contact results, the non-linear analysis gives useable 
results in the contact region. It was shown that reliable 
topographic information could be extracted from the con- 
tact data with sub-nanometer resolution. The data was 
most reliable at low drive velocities. 

A second material property obtainable in contact was 
the friction coefficient. The value obtained here at the 
lowest velocity, [i = 0.35, was consistent with values pre- 
viously obtained for Si-SiJ^ However, the results were 
sensitive to the point in contact at which it was mea- 
sured, or else to the range chosen if it was averaged in the 
contact region. Also, they deteriorated badly as the drive 
velocity was increased. It was concluded that the results 
were being effected by variable drag ) 12 i 13 which was ne- 
glected here. This effect is exacerbated in high viscosity 
liquids and at high drive velocities. Fortunately, it was 
shown that the pre-contact forces in the analyzed experi- 
mental data were not sensitive to the value of the friction 
coefficient. However, to improve the performance in high 
viscosity liquids of this method for measuring friction, 
variable drag will have to be accounted for. 
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